The Density Functional via Effective Action 
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A rigorous derivation of the density functional via the effective action in the Hohenberg-Kohn the- 
ory is outlined. Using the auxiliary field method, in which the electric coupling constant need not 
be small, we show that the loop expansion of the exchange-correlation functional can be reorganized 
so as to be expressed entirely in terms of the Kohn-Sham single-particle orbitals and energies. 

PACS numbers: 71.15.Mb 



Interactions among electrons largely determine the 
structure, phases, and stability of matter. Pragmatic 
advances in this subject, however, are nontrivial. When 
the number of electrons involved becomes large, calcu- 
lations based on constructing many-electron wave func- 
tions soon lose a,ccuracy and will be stopped by an "ex- 
ponential wall" 1] . Density functional theory (DFT) , us- 
ing the three-dimensional electronic density as the basic 
variable, is free from this wall. DFT originated from the 
theorem of Hohenberg and Kohn (HK)'2|, which states 
that there exists a unique description of a many-body 
system in its ground state in terms of the expectation 
value of the particle-density operator. The HK theorem 
assures that the ground state energy Eg is obtained by 
minimizing the energy functional E^j with respect to the 
electronic density n: 



Eg — min E^, [n] 



(1) 



Mermin [sj extended this theorem to finite-temperature. 

To make practical use of the HK theorem, a suitable 
computational scheme is necessary. Kohn and Sham [4] 
proposed a decomposition scheme, aiming to express 
E^j [n] via an auxiliary, noninteracting system that yields 
a particle density identical to that of the physical ground 
state. For a nonrclativistic fcrmion system described by 

H = I d^ij\^) ( -^V^ + w(x) - ^ j V;(x) 



^t(x)v3t(y)V3(y)Vi(x) 



dxdy, (2) 



the energy functional, with representing the electric 
coupling constant and To[n] being the kinetic energy of 
the auxiliary system, takes the form 



E,, [n] 



u(x) n (x) dx — /iiVe + Tq [n 
2 r !■ n(x)n(y) 



dxffy + E,j,^ [n] , (3) 



where \i = chemical potential, = number of electrons, 
v{k) = external potential, and E^c [n] is the so-called 
exchange-correlation energy functional. This exact de- 
composition cannot exist without the quantity ^-^""^"'^ 



being well defined. Being independent of w(x), the sum 
of the last three terms in ([3]) is universal. All of the many- 
particle complexity is now completely hidden in E^c [n] ■ 
Although TnWi] + Excin] admits no free parameter and 
is universal [2], its explicit construction remains elu- 
sive, and parameter-containing empirical functionals are 
therefore introduced. Cases of failure and limitations of 
these empirical functionals have been discussed [ill. On 
the other hand, a number of groups 0,BS[l3| have pur- 
sued first-principle derivation of the density functional 
via effective action. These efforts either introduce an 
auxiliary field 0, [13] or expand in powers of 0, 0]- 
The strengths of the auxiliary field approach are the sim- 
plicity of the effective action expression and the fact that 
each term already includes infinitely many Fcynman di- 
agrams fvi\. However, this approach seems ^] to lack a 
direct connection to the Kohn-Sham (KS) scheme. Such 
a connection can be made in the expansion in powers of 



[3, 12 , but that expansion is good only when 
small [13j . The validity of that assumption depends on 
the strength and variation of u(x). 

In this Letter, without assuming small, we report 
our development IJ] of an auxiliary field method that 
makes a direct connection to the KS scheme. To lighten 
the mathematical expressions in our finite-temperature 
formalism, we suppress the spin degree of freedom (as it 
is easy to include) and denote by a dot (circle) the three 
(four) dimensional integral contraction (with r denoting 
the Euclidean time, x = (t, x)) 



i-b = 



dx a(x) 6(x) 



aob = J dTdxa{T,x) b{T,x) — J dx a{x) b{x) . 

To probe the electron density, one introduces to H a 
classical source term J(x) coupled to '0t(x)7/'(x), H 
H + J-^ijj'^^fj) = Hj. Let (3 be the temperature inverse, 
l3J-{tfj'^tp) is written as J o^ip'^ ip) = J dxJ {x)'(j:i'^ {x)il}{x) . 
The partition function now is a functional of J, that is 



Z[J] 



-I3W[J]. 



Tr 



ETr 



(4) 



5n 



To disentangle the quartic fermionic interaction, we 
use the standard procedure of introducing an auxiliary 



2 



field (f> and express Z[J] as a path integral over both the 
Grassmann fields and the auxiliary field 



(5) 
where 

5[(/),V^,V'] = -^Ti-\n{u) + ^<j)ouo(t) + ^pKG-K'ilj (6) 
G-\x, x') = (9r + h{x) + i(uo<j>)^ + J(a;)) S{x - x') (7) 

/i(x) = hUion(x)-M (8) 

2m 

u{x, x') = (5(r - r')eV|x - x'| = (5(r - r')u(x, x') , (9) 

with -0'^^-' denoting the Grassmann fields satisfying 
tp^'^^il3,x) = — ■i/'^'l') (^0, x). It is easy to verify that 



SJ{x) 



{i;^x)ijix)}j^{nix))j=njix) . (10) 



Eq. pO)) expresses n in terms of J. The effective action 
is defined as the Legendre transformation of /3W^[J] 



T[nj] = PW[J] - J on J 



(11) 



where the subscript J indicates that the domain of T[n\ 
is the set of density profiles reachable by varying J. 
Eq. (fTT|) also leads to 

ST[n] 



-J 



(12) 



We now show that Ev[n] — lim^^oo ;gr[n]. Eq. ([4]) as- 
sures that at the zero temperature limit W[J] is simply 
the ground state energy corresponding to Hj. Evidently, 
when J = 0, lim/j^oo ;gr[n] |„=„^ = H^[J]|j=o = Eg 
where Eg stands for the ground state energy correspond- 
ing to H and Ug represents the electron density at the 
physical (J — 0) ground state. When J 7^ 0, the corre- 
sponding electronic density nj is different from Ug and 
lim^_^oo ^r[n]|„^„j represents the expectation value of 

H, calculated using the ground state wave function cor- 
responding to a different Hamiltonian Hj. Thus by the 
definition of the ground state, lim^^oo ^r[n] |„=„j > 
lim/3^oo ^r[n]|„=„g. This means that lim^^oo |r[n] 

reaches its minimum at Ug. Thus lim^^oo ^Ti*^] has all 
the properties attributed to the energy functional E^j in 
^ and ([3]). Since the HK theorem states that this func- 
tional is unique, it must in fact be equal to lim^^oo ^r["-]- 
If we make a change of variable — > -I- iu~^oJ in 
(l5][7l) and integrate over the Grassmann fields, we obtain 



-f3W[J] 



rjou Oj 



D(j)e 



where 



m = -iTrln(w) + ]^c^ouocj) - Trln(G^i) 



(13) 



(14) 



and 

G^^{x, x) = (dr + /i(x) + i(uo0)^) Six - x') . (15) 
Eq. implies that 

(16) 



(3W[J] ^ (3W4J] - -Jou-KJ, 



and thus the left-hand side of pO|) can be expressed dif- 
ferently, leading to 



nj 



iip ~ u J 



(17) 



where iifi = S{l3Wci,[J])/SJ . To evaluate f3Wc/,, we follow 
Jackiw [11] and let — > + (/j in (|13m5p . In particular, 
(fTSl) is rewritten as 



(18) 



G^l^{x, x') = G-^{x,x') + iS{x - x') (uocj))^ , 
and one obtains 11] 

PW^J] = iTrln(p-iou) + ^^ouo^ - Trln (G-^) 



00 ^ 

-iJo(p — > — 



^/(''^'Mo6i...o6fc 



/IPI, conn 



,(19) 



where the subscript "IPI, conn." means to include only 
connected, one-particle-irreducible diagrams, b = uocj), 



f)-^ = u-^ - D 



D{x,y) = G^{x,y)G^{y,x) , 



and 



_ (-1) 



fc-i 



dxi . . . dxk 



G^{xk,xi) . . .G^{xk-i,Xk){ib{xi)) . . .{ib{xk)) ■ (20) 

Fukuda et al. [3] obtained an expression similar to (fT9|) 
and used it to derive an effective action as a functional 
of if. They also noted that this auxiliary field approach 
does not make a direct connection to the KS scheme. 

Coming to the point of departure from typical auxiliary 
field approaches, we show below how an exact correspon- 
dence to the KS scheme can be made for the auxiliary 
field method by decomposing the source J in a particu- 
lar way. Let us define a free fermion propagator Qo by 

go\x,x') ^ \dr + /i(x) + Jo{x)] S{x - x') , (21) 

where Jo is chosen (if ^%^^lnj exists, Jq exists and can 
be written as u-nj + ^%^^|tij + J ) such that 



Sn I"./ ' 
Goix^x) = nj{x] 



(22) 



Eq. demands that this non-interacting (KS) system 
have electron density, —Goix,x), identical to nj{x), the 
electronic density of the physical system (where Coulomb 



3 



interactions exist). In (jl9p . each occurrence of iuoip 
through Gtp is to be replaced by J + uonj (from HT])). 

To bring out the KS scheme, we perform the following 
source decomposition 

J[n] = {Jo[n] ~ uonj) + J'[n] = Jo[n] + J'[n] . (23) 

Then from (fT5|) and ([TT]) we have 

G-\x, x') = g^^ix, x') + J'{x)S{x - x') . (24) 



Although the source decomposition (|23p is introduced 
here for the first time in the auxiliary field approach, a 
similar method was used in 0] to perform perturbative 
calculations using as the expansion parameter. 

Substituting ((T7]) and (dH) into one obtains an 

expression for /3M^[J], which, upon introducing a param- 
eter A (to be set = 1 in the end) to denote the loop order, 
has the form f3W[J] =JWo[J] + J2Zi y{PW,[J+uonj]), 
where in particular [l4| 



(3Wa[J] = (3Wo[J + Uonj] - -njououj 



(25) 



with pWo[J + uonj] = -Trln(G-i). 

To arrive at an expansion headed by — Trln(^Q'^) 
instead of —Trln{G^^), and containing the expression 
Wi [Jo] instead of Wi [J+uonj] , we expand Wi [J+uonj] — 
Wi[Jo + J'] in powers of J' (subscript I omitted in the 
equation below) 



W = W[Jo] 



SWjJo] 
SJq 



1 S^WjJo] 

2 6Jo SJq 



ij'oj' 



(26) 



The expression Wi [Jq] means that J is replaced by Jq but 
Uonj is kept unchanged. 1J| With l|26p. we may express 
/3VF[J] as a double series 

I3W[J] = (3Woo +(3J2w^kil- <5«.o4,o) J''' A' , (27) 

i,k 

where each Wik involves the fc'th derivative of Wi. In 
particular, Wqo is given by (with nj n hereafter) 



pWoo = pWoo 
and in view of 



1 



-Trln(g, 



Wqi is given by 



1 



nouon , (28) 



S{f3Wo[Jo]) S{f3WoQ[Jo] 

= n = 



SJo SJo 
The second half of (|29|) suggests that we define 
fo[n] ^ pWooiJo] - Joon , 



(29) 



(30) 



the Legendre transformation of the zeroth order contri- 
bution from /3VK[J] (in terms of J' and A), leading to 



Sn 



-Jo 



(31) 



Comparing ([3T|l with p2|) . we find 



5{r[n]-To[n]) 



= -J' . 



(32) 



The idea now is to develop a series for r[n] led by fo[n]. 
Subtracting ([50)1 from pTjl . we have 



T[n] ~ foN = (3W[J] - pWooiJo] - J' on , 



(33) 



in which the last two terms on the right hand side exactly 
cancel the terms in Woo and Woi contributing to /3W[J]. 
So the series for F — Fq is just P7)) with those two terms 
removed. Next we convert the double sum in (j27p into a 
single sum by expanding J' as a series in A. We write 



oo 



J'[n] ^^Ji[n]X 
1=1 



(34) 



where the precise expressions for Ji , J2 , ■ ■ • are as yet 
undetermined since p4p is not a loop expansion. We 
substitute (p4)) formally into ((33)) and (p7)) to obtain a 
series 



F[n]-fo[n] 



ran] A' 



(35) 



in which each F; is defined explicitly in terms of the J^, 
PWk<i[JQ]j and their derivatives. Because Wqi is missing 
from ([55]) . any occurrence of Jk is accompanied by at 
least one other factor Jk' or else by an occurrence of some 
Wi>o, and hence by a power of A higher than the /c'th. 
In other words, the expression for F/>i involves only Jj- 
with k < I. We finally remove the indeterminacy in ((34|) 
by imposing ()32p to hold order by order in A, leading to 



STi[n] 
6n 



= -Ji 



(36) 



Since F/>i involves only Jk<i, all the J; and F/ 
can be found explicitly by applying ([55]) and 
alternately. The first few expressions are Fi 
pW^[Jo] = -'^TrHV-^j^ou), Ji = -^^^S^< 



■5 Jo 
Sn 



PW2[Jo] 



SJo 



■Jl 



.Jl. 



1 T S'{PWo[Jo]) 
2"'1° SJoSJo 

For an arbitrary Jq, one will obtain a corresponding 
density n. The computation of ^F[n] using (|5D|) . and 
(|36p evaluates the energy functional at density n, which 
may or may not be the ground state density. To obtain 
the ground state density and the corresponding Jq, one 
needs to solve at zero temperature limit the extremal 



equation 



sr[n ] 

Sn 



, which we turn to shortly. 



To carry out the calculation of F[n], we need to com- 
pute Jl (see ([36)) ) via the functional derivative 



S_ 

Sn 



Sn 
SJ^ 



SJn 



SJo 



(37) 



Diagrams corresponding to /3W/[Jo] and their derivatives 
contain the u, Qoj and Do = Vj^j^^ propagators. It is 



4 



easy to show that one may express Sn{x)/5Jo{y) as 



5 My) 



Qoix, y) Goiy, x) ^ Dj^j (x,y) (38) 



and thus ^ ~ D = , which we cah the inverse density 

correlator. The differentiation rules of Qq, Vq, and Dq"^ 
with respect to Jo can be expressed diagrammatically: 



x' ^ H. X u(x, x') 

x' . ^ 



X D^\x,x') 



X Qq{x,x') x' tANVVVW* X Vq(x,x') 

XX 



6go{x, x') _ d 
6Jo{y) SJoiy) 



S'Do{x,x') _ 6 
SJo{y) SJoiy) 



= - \y 

x' 

X 



X 
X 



SDo\x,x') 
SJoiy) 



SJoiy) 



The differentiation rules of Qo, T^o, Siiid Dq^ with re- 
spect to n are simply obtained by compounding the re- 
sults above with (IST)) . We show only one example: 



Sgoix,x') S 



Sn^z) Sniz) 



Equipped with these differentiation rules, one may 
use standard diagrammatic expansion to compute the 
Wi [Jo]s, their functional derivatives with respect to Jo, as 
well as JiS to facilitate the calculations of Tis. Because 
Doix,y) = go{x,y)goiy,x), both Vq = (u^^ ~ Dq) 
and can be expressed in terms of single-particle or- 
bitals and energies through Goix^y) -the propagator of 
the KS system- which can be expressed as 



Goix,y) = ^0a(x)(?!);(y) 



g (ea-M)(Tx— Tj,) ^ 



i~na) if Tj. < Ty 



(39) 



(1 - Ha) if Tx > Ty ' 

where tIq = l/(e'''^'^"^^^-|-l), ^^Ua = N^, and the single 
particle orbital 0q(x) satisfies 

h{x.) + Jaix) (l)a{x.) ^ (ea ~ fJ.)(f>ai^) ■ 



Since 



= — Jo, the extremal condition = ^^^-M 



that determines Ug and Joi^^g] (as /? —^ oo) becomes 



Sn 



'5(E,=ir.M)_ sij:z,r,[Jo[n] 



Sn 



SJo 



Jo- (40) 



Eq. (j40|l has to be solved self-consistently by keeping Ti 
terms up to some order in A. Although a truncation is 
necessary, we note that each diagram in our expression al- 
ready corresponds to infinitely many Feynman diagrams 
when using as the expansion parameter. This is easily 
seen by performing the small expansion of Vq 

Vq — U + UoDqoU + UoDqoUoDqoU -|- . . . , 

a sum of infinitely many (dressed) propagators. Interest- 
ingly, in the strong coupling limit where one must treat 
as a small parameter, we may express I?o sls 



-Do 



-Do'- 



°^o ' 



while the traditional expansion fails completely. 

Finally we sketch how ^ arises from V[n\. Eq. ((30|) 
may be rewritten as -gFoM = ^ [^Tr ln(t?Q""'") — Joon] -f- 
j^nouon. Because —Tiln{QQ^) = '^^^^^{1 — Ua), at 
zero temperature limit, the first two terms of ^Fq 
above give rise to the To[n] — iiN^ + J u(x)n(x)dx 
while the last part is exactly the Hartree term [14| . 
The exchange-correlation functional Exc[n] equals 
lim/3_oo i Fi[n]. We also comment that the ex- 



citations of the system can be studied under this 
formalism and the energy functional shown in this letter 
has the correct single-electron limit 1J| . 

Providing a scheme beyond perturbative expansion in 
, we have proposed an effective action construction that 
will contribute to the development of the parameter-free 
universal density functional. 

This research was supported by the Intramural Re- 
search Program of the National Library of Medicine of 
the National Institutes of Health. 
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